Friction effects and clogging in a cellular automaton model for pedestrian dynamics 
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We investigate the role of conflicts in pedestrian traffic, i.e. situations where two or more people 
try to enter the same space. Therefore a recently introduced cellular automaton model for pedestrian 
dynamics is extended by a friction parameter ^. This parameter controls the probability that the 
movement of all particles involved in a conflict is denied at one time step. It is shown that these 
conflicts are not an undesirable artefact of the parallel update scheme, but are important for a correct 
description of the dynamics. The friction parameter ^ can be interpreted as a kind of internal local 
pressure between the pedestrians which becomes important in regions of high density, ocurring e.g. 
in panic situations. We present simulations of the evacuation of a large room with one door. It is 
found that friction has not only quantitative effects, but can also lead to qualitative changes, e.g. 
of the dependence of the evacuation time on the system parameters. We also observe similarities to 
the flow of granular materials, e.g. arching effects. 

PACS numbers: 45.70.Vn, 02.50.Ey, 05.40.-a 
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I. INTRODUCTION 

Methods of physics and modern computer science have 
been used successfully for the investigation of vehicu- 
lar traffic problems for a long time ^, |^. Also 
pedestrian dynamics has attracted some attention in re- 
cent years Q and many interesting collective effects and 
self-organisation phenomena have been observed (for an 
overview see ||, ^, ^). 

One topic which has not been studied intensively up 
to now (see, however, is the relevance of local con- 
flicts in pedestrian traffic. A conflict indicates a sit- 
uation in which two or more people try to enter the 
same space in one timestep. Obviously this is a real 
two-dimensional effect which has no counterpart in (di- 
rected) one-dimensional vehicular traffic. These conflicts 
are local phenomena which can have a strong influence on 
global quantities like evacuation times and flows in the 
presence of bottlenecks. Typical examples where con- 
flicts become important are situations with clogging and 
stucking encountered in crowds of panicing pedestrians, 
e.g. near intersections and bottlenecks. In real life this 
often leads to very dangerous situations and injuries or 
even fatalities during evacuations. From a physics point 
of view these phenomena are related to flow properties 
of granular materials 0, ^, ^, |l^ . 

Depending on the choice of parameters, the cellular 
automaton (CA) model introduced in |5[ |l6j is able 
to reproduce many of the observed collective effects of 
pedestrian dynamics, i.e. panicing and herding behaviour 
[ p^ . In this paper the model is extended by a friction 
parameter /.t, which allows an improved description of 
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clogging and stucking phenomena of pedestrians. It will 
be shown that conflicts and friction are responsible for 
several interesting effects which become especially rele- 
vant in high density situations, e.g. during evacuation 
processes, and that they are important for a correct re- 
production of the dynamics. We start with a short sum- 
mary of the model's basic concepts and its update rules. 



II. BASIC PRINCIPLES OF THE MODEL 

The model considered here is a CA where the space is 
discretized into small cells which can either be empty or 
occupied by exactly one pedestrian. Each of these pedes- 
trians can move to one of its unoccupied neighbour cells 
(see flg. 0) at each discrete time step t ^ t + 1 accord- 
ing to certain transition probabilities. The probabilities 







p 

-1,0 





p 

0,-1 


P 

0,0 


P 

0,1 





p 

1,0 






FIG. 1: Deflnition of the transition probabilities pij. 



are given by the interaction with two discrete floor fields^ 

and 
and 



D and S P, P, |T|, 0. The field strengths D,^ 
Sij at site (ijj) are interpreted as number of D- 
S-particles, respectively, present at that site. The two 
fields determine the transition probability in such a way 
that a particle movement in the direction of higher fields 
becomes more likely. 

The dynamic floor field D represents a virtual trace 
left by moving pedestrians. Similar to the process of 
chemotaxis ^ used by some insects for communica- 
tion, this trace has its own dynamics, namely diffusion 
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and decay. Two parameters, a £ [0,1] and S € [0,1], 
control the broadening and dilution of the trace. Every 
moving pedestrian creates a D-particle at its origin cell. 
Since there is no restriction of the maximal number of 
D-particles at a site, D can be regarded as a bosonic 
field. 

The static floor field S, on the other hand, does not 
change in time. It reflects the surrounding geometry and 
e.g. specifying attractive space regions. In the case of 
the evacuation processes considered here, the static floor 
field describes the shortest distance to an exit door. S is 
calculated for each lattice site using some distance metric. 
The field value increases in the direction of the exit such 
that it is largest for door cells. An explicit construction 
of S can be found in [|7[|g|. 

The floor fields are used to translate a long-ranged spa- 
tial interaction into a local interaction, but with memory, 
similar to the phenomenon of chemotaxis in biology. The 
only other model which so far reproduces all observed col- 
lective effects of pedestrian flow, the social-force model 
pl| , uses exponentially decaying repulsive forces between 
pedestrians. In contrast, in our approach the interac- 
tion is local and attractive. However, pedestrians do not 
interact with the density, but with the velocity density 
created by the other particles. 



A. Update rules 

The update rules of the full model, including the inter- 
action with the two floor fields, have the following struc- 
ture: 

1. The dynamic floor field D is modified according 
to its diffusion and decay rules [0, controlled by 
the parameters a and S. In each time step of the 
simulation each single boson of the whole dynamic 
field D decays with the probability 6 and diffuses 
with the probability a to one of its neighbouring 
cells. 

2. For each pedestrian, the transition probabilities pij 
for a move to an unoccupied neighbour cell 

(fig. |l|) are determined by the local dynamics and 
the two floor fields. The values of the fields D and S 
are weighted with two sensitivity parameters ks G 
[0,oo[ and kjj € [0, oo[. This yields 

Pij = exp {koDij) exp (fcs5y )(l - , (1) 

with the occupation number = 0, 1, the obstacle 
number 



for forbidden cells (e.g. walls) 
else 



and the normalization 

A^ 



(2) 



(3) 



3. Each pedestrian chooses randomly a target cell 
based on the transition probabilities pij determined 

by (0). 

4. Conflicts arising by any two or more pedestrians 
attempting to move to the same target cell are re- 
solved by a probabilistic method. The pedestrians 
which are allowed to move execute their step. The 
explicit procedure of conflict resolution will be de- 
scribed in detail in section III . 



5. I? at the origin cell (i, j) of each moving particle is 
increased by one: Dij Dij + 1. 

The above rules are applied to all pedestrians at the same 
time (parallel update). This introduces a timescale of 
about 0.3 sec/timestep [|4|. 



III. RESOLUTION OF CONFLICTS 

Due to the use of parallel dynamics it is possible that 
two or more particles choose the same destination cell in 
Step 3 of the update procedure. Such situations will be 
called conflicts. At first it appears that conflicts are un- 
desirable effects which reduce the efficiency of simulations 
and should therefore be avoided by choosing a different 
update scheme. We will show in the following that this is 
not the case and that confiicts are important for a correct 
description of the physics of crowd dynamics. 



Conflicts without friction 



In [|14|, |17| conflicts between pedestrians were solved in 
the following way: whenever m > 1 particles share the 
same target cell, one (/ € {1, . . . , m}) is chosen to move 
while its rivals for the same target keep their position. 
There are two main ways to decide which particle / is 
allowed to move: 

1. According to the relative probabilities with which 
each particle chooses its target cell, i.e. the proba- 
bility for particle I to move is pfj/J2"=iPif- 

2. All particles move with the same probability ^. 

The observed behaviour has been shown to be quite ro- 
bust and does not depend on the details of the conflict 
resolution [H. 



B. Introduction of the friction parameter /i 

We now extend the basic model by a new friction pa- 
rameter /i e [0,1], in order to describe clogging and 
stucking effects between the pedestrians. Whenever two 
or more pedestrians try to attempt to move to the same 
target cell, the movement of all involved particles is de- 
nied with the probability /i, i.e. all pedestrians remain 
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at their site (see fig. g). This means that with proba- 
bility 1 — fJ- one of the individuals moves to the desired 
cell. Which particle actually moves is then determined 
by the rules for the resolution of conflicts as described 
in Sec. [II A, Note that for vanishing friction /i 
recover the dynamics studied in 0, ^ |l7| . 



we 




t+ 1 



FIG. 2: Refused movement (for m 
parameter ^. 



4) due to the friction 



The friction parameter might be interpreted as the ef- 
fect of a moment of hesitation: Pedestrians in conflict sit- 
uations slow down or hesitate for a short moment when 
trying to resolve the conflict. This reduces on average 
the velocities of all involved particles. 

With the above definition of fi and the extended up- 
date rule it is easy to see that fi works as some kind of 
local pressure between the pedestrians. If fi is high, the 
pedestrians handicap each other trying to reach their de- 
sired target sites. As we will see, this local effect can 
have enormous influence on macroscopic quantities like 
flow and evacuation time. We like to point out that fric- 
tion as introduced here does not reduce the velocity of a 
freely moving particle. Its effects only show up in local 
interactions. 

As one can see it is necessary for such kind of inves- 
tigations to use a parallel update in the model. Any 
other form of random or ordered sequential update will 
disguise the real number of arising conflicts between the 
pedestrians in the system. 



IV. SIMULATIONS AND RESULTS 

In the following we describe results of simulations for 
a typical situation, i.e. the evacuation of a large room 
(e.g. in case of fire). In it is shown that by variation 
of the sensitivity parameters ks and ko (see Sec. [I A) 
three main regimes for the behaviour of the particles can 
be distinguished. For strong coupling ks and very small 
coupling ko we find an ordered regime where particles 
only react to the static floor field. The behaviour then 
is in some sense deterministic. The disordered regime 
characterized by strong coupling k]j and weak coupling 
ks leads to a maximal value of the evacuation time T. 
Here the behaviour is typical for panic situations. Be- 
tween these two regimes an cooperative regime exists 
where the combination of interactions with the static and 
the dynamic floor fields minimizes the evacuation time. 



Fig. g shows a schematical phase diagram in the space of 
the couplings ks and ku. 



disordered 




ordered 



FIG. 3: Schematic phase diagram. 

Here we want to focus on the influence of the friction 
parameter /i on the evacuation times in the three main 
regimes. As we will see the influence of /x is strongest in 
the ordered regime. This will lead to a new interpretation 
of this regime. 

We consider a grid of size 63 x 63 sites with a small exit 
of one cell in the middle of one wall. The particles are 
initially distributed randomly and try to leave the room. 
The only information they get is through the floor fields. 
The strength of the static floor field at a site («, j) is in- 
versely proportional to the distance of (i, j) from the exit 
measured using the metric described in [p^ . Fig. ^ shows 
a typical stage of the dynamics for an initial particle den- 
sity of p = 0.3, i.e. 1116 particles. A typical half-circle 
jamming configuration in front of the door develops. 




FIG. 
tion. 



4: Typical middle stage of the dynamics of an evacua- 
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A. Influence of /i and p 

First we look at the averaged evacuation times T in the 
three main regimes, in dependence of the particle density 
p and the friction parameter /z. All evacuation times are 
averaged over 500 samples and measured in update steps. 
Figs. ^ and ^ show the influence of a varying /i parameter 
on the three regimes for the low density p = 0.03 and the 
high density p = 0.3. 



— ordered 

— — disordered 

— cooperative 




0.0 0.2 0.4 0.6 



way to the exit. Even for low densities they will form 
a jam in front of the door after a short time. But since 
these jams are small, the evacuation time increases only 
for large friction parameters p> ^. 

The behaviour is different in the high density regime 
p — 0.3 (fig. H). One finds a weak increase of the evacu- 
ation time in the disordered and cooperative regime for 
high p values (p > 0.6). Even if the particles are not 
packed close together in front of the door, they form a cue 
and hinder each other. Increasing p leads to a sharp in- 
crease of the evacuation time in the ordered regime which 
first becomes larger than that in the cooperative regime 
and finally exceeds even that in the disordered regime. In 
fact it diverges for p — > 1. This behaviour can be un- 
derstood from the microscopic configurations occuring. 
A short time after the start of the evacuation nearly all 
particles of the system are forming a big jam in front 
of the door. In this large density region many confiicts 
occur and for large values of p the outfiow is strongly 
surpressed. Here the pressure between the pedestrians 
becomes so strong that any motion is almost impossible. 
Such a behaviour has been observed in panic situations 
and also in simulations using the social-f orce m odel 



It will be discussed in more detail in Sec. |VB. 

These results are supported by figs. ^ and ^ which show 
the infiuence of an increased density p for fixed values of p 
on the evacuation times. In the disordered regime even 



FIG. 5: Average evacuation times T for a large room in depen- 
dence of the friction parameter p for the low density p = 0.03. 
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FIG. 6: Same as fig. but for a higher density of p = 0.3. 

In the low density regime, p = 0.03 (fig. ||), increasing 
p has only a very weak effect on the evacuation times 
of the disordered and the cooperative regime. Here the 
pedestrians move almost independently of each other so 
that confiict situations, even close to the door, are rare 
since almost no jamming occurs. The behaviour is differ- 
ent in the ordered regime. All particles find the shortest 



FIG. 7: Density-dependence of the evacuation time in the 
ordered and disordered regime for p = 0.9. 

the infiuence of the large value p — 0.9 is only visible 
for higher densities (p > 0.1). The cooperative regime 
shows a similar behaviour. In contrast, in the ordered 
regime even for very small densities p the influence of 
the high value p — 0.9 is very strong (fig. |^). This is 
supported by fig. ^, where the relative time differences 
for the two values /i = and p — 0.9 for all three regimes 
are shown. Here the very strong increase of the relative 
time difference in the ordered regime after a very small 
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FIG. 8: Relative time difference of the evacuation times for 
/i = and fi = 0.9 in the three regimes. 



FIG. 9: Evacuation time as function of the sensitivity param- 
eter ks for different fi values (ko = 0, p — 0.3). 



density of p = 0.003 is remarkable. After this increase 
the relative time difference remains almost constant. 



B. Reinterpretation of the ordered phase 

The introduction of the friction parameter /i allows an 
improved interpretation of the ordered regime, i.e. the 
case of strong coupling {ks > 3) to the static field S and 
weak coupling to dynamic field D. This regime is almost 
deterministic, with pedestrians moving straight towards 
the exit on the shortest path. Now an increased ^ value 
introduces a negative interaction between the particles 
into the system, i.e. the pedestrians hinder each other 
due to strong competition for the unoccupied target sites 
near the exit. So a strong coupling to S together with a 
high fi value, which works as an internal local pressure 
between the particles, describes a typical panic situation, 
where an ordered outflow is inhibited due to local con- 
flicts near bottlenecks or doors, resulting in strongly in- 
creased evacuation times (such situations are well known 
from emergency evacuations due to fire or other reasons 
in sports arenas or passenger vessels). This can be seen 
in fig. ^, where the influence of an increased coupling 
strength to S for fixed /i is shown. For ks — > the par- 
ticles perform a pure random walk and the evacuation 
times are very large and almost independent of /z. In 
this situation conflicts between the particles are not very 
important for the dynamics. In contrast, for ks — > oo 
they choose the shortest way to the exit. Since for /i = 
there is no internal pressure between the particles the 
evacuation time is minimal. However, for /i — > 1 the 
number of unsolved conflicts increases with ks due to 
the strong jamming at the exit. This results in stucking 
and clogging phenomena and highly increased evacuation 
times. For very high fi values (/i = 0.9 in fig. O) one finds 



a minimal evacuation time for an intermediate coupling 
(ks ~ 1). This means that a larger ks, which implies a 
larger average velocity of freely moving pedestrians, leads 
to larger evacuation times. This collective phenomenon 
is very similar to the faster-is-slower effect ||, ||]. Note 
that a similar minimum is characteristic for the coop- 
erative regime. However, there local minima have been 
observed as a function of the coupling ko to the dynamic 
floor instead of ks [ p7[ . 



C. Time evolution of an evacuation 

In the following we look at the time evolution of an 
evacuation, i.e. the number of people who left the room 
at a certain time stage. Figs. |lo[and show the time de- 
pendence of the number of evacuated persons for all three 
regimes and different fi values for a high density p = 0.3. 
The curves show a nearly linear increase since the very 
high initial density leads to strong clustering at the door 
already at the beginning of the evacuation. The evacua- 
tion times are strongly increased due to the large value 
/i = 0.9. This effect is again strongest in the ordered 
phase. For fi — the corresponding evacuation time is 
the smallest whereas for /i = 0.9 it becomes the largest 
one of all regimes. At the end of the evacuations, when 
most of the particles have left the room, the gradient of 
the curves becomes rather small. 

Because of the averaging over many samples the curves 
are very smooth. To have an impression of the evolution 
of one single evacuation and of the variance of evacuation 
times, fig. ^ shows the averaged curves envelloped by 
the curves with the minimal and the maximal evacuation 
time of the whole sample for the ordered phase and three 
fi values. For fi = the evacuation process is nearly 
deterministic in the ordered regime and the fluctuations 
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FIG. 11: Evacuated persons in dependence of time for all 
three regimes; p = 0.3 and p — 0.9. 

(due to the random initial conditions and the dynamics) 
are very small. With increasing /x values the internal 
pressure is increased and the envelloping curves differ 
clearly from the averaged curves. 

In figs. |l^ and ^ again the averaged curves together 
with the extremals are shown for the two densities p = 
0.03 (fig. |l|) and p = 0.003 (fig. |l|) and the two friction 
parameters /i = and p, — 0.9. Here the same effects as 
for p = 0.3 can be seen, but much clearer. 

The time evolution of one single sample exhibits an 
interesting dynamics. In figs. |l2|-|l^ small plateaus can be 
observed where over short time periods no persons leave 
the room. This irregular behaviour is well-known from 
granular flow and is typical for clogging situations ||^, 

|l^ ^ |l^ . The plateaus are formed stochastically 
and can therefore not be observed after averaging over 
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FIG. 13: Evacuated persons in dependence of time for the 
ordered regime; averaged, longest and shortest evacuation 
process for two p values {p = and p = 0.9) and density 
p = 0.03. 



various samples (figs. ^ 

Note that also the variance of the evacuation time in- 
creases strongly with p. Fig. ^ shows that the average 
evacuation time T is not always a meaningful quantity 
for safety estimates since the variance can become quite 
large. For the very small density p = 0.003 the gradient 
of the curve is rather fiat at the beginning of the evacu- 
ation (fig. |lj). Here there are only few particles in the 
system and no cue is formed near the exit. 
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FIG. 16: An example for the time evolution of the number of 
evacuated persons. 



FIG. 14: Evacuated persons in dependence of time for the or- 
dered regime; averaged, longest and shortest evacuation pro- 
cess for two values {fj, = 0, fj, = 0.9) and density p = 0.003. 



D. 



Mean-field approximation for the ordered 
regime 



As we have seen in Sec. IV C the curve of the num- 



ber of evacuated persons N{t) grows almost linearly in 
the ordered regime, especially for high densities p. For 
this regime we now calculate approximatively the /z- 
dependence of this curve, i.e. N = N{t,fi). 

As explained earlier, in the ordered regime after a short 
time a big jam forms at the door due to the strong cou- 
pling to the static field S. Fig. ^ shows a typical local 
configuration in front of the exit. Generically, 3-particle 



Door 



N(t) 



FIG. 15: Typical configuration in front of the door in the 
ordered regime. 

conflicts over the unoccupied lattice site in front of the 
door occur. At time t, with probability 1 — /i one of these 
particles is able to move. In the next time step t + 1 this 
particle will escape through the door with probability 1. 
Neglecting conflicts elsewhere in the system, because of 
the big jam in front of the door the configuration shown 
in fig. |l5| will be restored at time t + 2. Therefore, re- 
peating the above sequence a typical representation of 
the time evolution of N{t) is shown in fig. ^ 



Using simple combinatorics, the expectation value of 
N{t) can be written as 



lt/2] 



c 

where C is the normalization factor defined by 

[t/2] 



(4) 



t - n 



(1 - 



[t/2] 



1- (a^- l)*+i 
2 - ^ ' 



(5) 



In the last step we have used the identity (Aq) which is 
derived in Appendix Due to (^) and (d|7 (^(i)) is 
related to C by the equation 

'^-C= - C+^^:^^ {N{t)) C (6) 



or 



d/i /i /i(l — /i) 



(7) 



Finally an analytical expression for {N{t)) is obtained as 



{N{t)) 



(^-2)2 



0, (m-2)(/x-1)V(1 + ^) \ 
Ai(l + t)~2t ) ■ 

Asymptotically, (H) implies for large times 



{N{t)) - ^ — J- -t as t ^ oo 



(9) 



i.e. a linear behaviour as observed in simulations for the 
ordered regime. This expression is also consistent with 
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FIG. 17: Evacuated persons in dependence of time for the 
ordered regime; comparison of the numerical results with the 
analytical approximation IM) for the density p — 0.3. 
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FIG. 18: Time-dependence of the number of evacuated per- 
sons for the ordered regime; comparison of the numerical re- 
sults with the analytical approximation for the density 
p = 0.003. 



the fact that T ^ cx) for /i ^ 1. From fig. we see 
that the analytical approximation (^) for {N{t)) is very 
good for high densities {p = 0.3) and friction parame- 
ter values /i < 0.6. However, for very high parameter 
values (/X > 0.8) the agreement is not satsifactory. Here 
the assumption that a typical local configuration looks 
like that in fig. |l5| no longer holds. The observed flow 
is underestimated since the large friction parameter /i is 
responsible for a flow reduction through conflicts occur- 
ing away from the door. Therefore particles are not so 
dense-packed at the exit (compared to the situation at 
smaller /i- values) which leads to an effective reduction of 
conflicts. 

In the low density regime p — 0.03 (not shown in the 
figure) the approximation again gives satisfactory agree- 
ment for /I < 0.6. Only in the region of very low densities 



(e.g. p = 0.003, see fig. 18) strong deviations can be ob- 
served for all /i values. In this case there are only very 
few particles in the system, not hindering each other. 
In most cases they reach the exit independently and the 
evacuation time essentially depends on the diameter of 
the room (or the average walking length) rather than on 
the number of particles in the system. 



E. Number of conflicts during evacuation 



1 : disordered 
2: cooperative 
3: ordered 




FIG. 19: Number of conflicts in dependence of time t for all 
three regimes (p = 0.3 and p — 0). 



1 : disordered 
2: cooperative 
3: ordered 




FIG. 20: Number of conflicts in dependence of time for all 
three regimes (p = 0.3, p — 0.9). 

Next we investigate the time evolution of the number of 
conflicts arising during an evacuation process. Figs. |l^ 
and ^ show the number of conflicts in dependence of 
time for all three regimes in the high density case p = 0.3. 
The maximal number of conflicts occurs for all regimes 
and p values shortly after the beginning of the evacu- 
ation. The number of conflicts is always largest in the 
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disordered regime. As argued earlier here the motion 
resembles a random walk. Since there are many parti- 
cles in the system most of them will compete for empty 
sites with other particles crossing their path anywhere in 
the room. In the ordered regime all pedestrians gather 
around the exit quickly with most of them standing in 
the cue without competing [ p3| . Even though the num- 
ber of conflicts is smallest in the ordered regime, the in- 
fluence of the friction parameter on the evacuation time 
is strongest. The reason behind this is that nearly all 
conflicts take place in front of the door. Therefore they 
have a direct effect on the outflow which is reduced con- 
siderably. Conflicts arising away from the exit do not 
have a direct influence on the evacuation time. In fact 
they can even lead to an increased outflow b ecaus e they 
might surpress clogging at the door (see Sec. IVF). 








a very strong increase of the conflict density in the or- 
dered regime at the end of the evacuation for fj. — 0. 
During most of the time of the evacuation the majority 
of particles are forming a big jam and do not contribute 
to the conflict density. At the end of the evacuation the 
jam nearly is dissolved and only few particles are left in 
the room near the exit conflicting about the unoccupied 
sites. 



■ 1 : disordered 

■ 2: cooperative 




FIG. 23: Number of conflicts in dependence of the time for 
all three regimes; p = 0.03 and /i = 0. 



FIG. 21: Conflict density in dependence of time t for all three 
regimes (p — 0.3 and fj. = 0). 




FIG. 22: Same as fig. but with p = 0.3 and = 0.9. 

In figs. ^ and ^ the corresponding conflict densities, 
i.e. the number of conflicts divided by the number of par- 
ticles inside the room, are shown. In flg. Ell one can see 
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FIG. 24: Same as fig. but for p = 0.03 and p = 0.9. 

Figs. ^ and ^ show the number of conflicts and 
the conflict density during an evacuation for the density 
p = 0.03 and p, = 0.9. Since there are now less particles 
in the system, they will not meet each other that often in 
the disordered and cooperative regime. For times smaller 
than the evacuation time of the ordered regime the num- 
ber of conflicts and the conflict density are largest in the 
ordered regime, where all particles are packed in a small 
region near the exit. 
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F. Column in front of the exit 

As an example for safety estimations in architectural 
planning we investigate how evacuation times change if 
a column (i.e. a non-traversable obstacle) is placed in 
front of the exit. The authors of |j, who studied the 
same situation and found surprising results, called for 
complementary data from experiments or other models to 
confirm their findings. We therefore study the scenarios 
of 1^, 0]_in the following. 

Fig7^25| shows a mid time stage of evacuations for a 
column placed central in front of the exit. The size of 
the column is 3 x 3 cells. It is placed within a distance 
of one cell from the door. We compare with situations 



■ no column 
central column 

- column, shift one site 

■ column, shift two sites 





FIG. 25: Mid time configuration of an evacuation of a room 
with one door at the middle of the top wall with an additional 
column of size 3x3 cells placed in front of it. The width of 
the door is one cell. 

where this column has been shifted (to the left or right) 
parallel to the exit by one or two lattice sites. Our fo- 
cus is on the ordered regime where friction effects are 
strongest. The static floor fields S used in the simula- 
tions have been calculated using a Manhattan metric (for 
details see ) . Fig. ^ shows the averaged evacuation 
times as a function of the friction parameter fi for a room 
with different configurations of the column. The placing 
of the column docs not influence the evacuation times for 
/X = 0. Here the movement is nearly deterministic and 
the particles do not hinder each other. For higher val- 
ues of fj, the column becomes more and more relevant. It 
causes a subdivision of the crowd and decreases the local 
pressure between the particles. Therefore the evacuation 
times become smaller. They are minimized for a column 
placed in a slightly asymmetric configuration in front of 
the exit (with a shift of one lattice site, fig. p6|) . 

How can this surprising result be understood? On 
one hand, the column has a certain screening effect that 
forces some pedestrians to take a detour and therefore 
potentially increases the evacuation time. On the other 
hand, the column subdivides the pedestrian flow and so 
can lead to a reduction of conflict situations, especially 
close to the exit. The competition between these two ef- 
fects is then responsible for the nontrivial dependence of 
the evacuation times on the position of the column. 



FIG. 26: Evacuation time as function of the friction param- 
eter /i for four room geometries: no column, central column 
and column shifted by one and two lattice sites. The initial 
density is p = 0.3 and the coupling strengths are ks = 10, 
ko = corresponding to the ordered regime. 




FIG. 27: Number of conflicts for a room with one door com- 
pared to various positions of an additional column close to the 
exit. The friction parameter is fi = 0.9. All other parameters 
are the same as in fig. 



Fig. 27 shows the number of conflicts for three different 
positions of the column in comparision to a room with- 
out column. The evacuation time T is minimal with the 
column shifted one site (curve 3). In this case, the total 
number of conflicts during evacuation, which is the inte- 
gral of the curve, is also minimal. However, in the early 
stage of the evacuation, the number of conflicts is bigger 
than for all other cases. This is a result of the subdivi- 
sion of people in the early stage of the dynamics which 
occurs most pronounced in case 3. This subdivision en- 
ables pedestrians to escape from the room faster than in 
the other cases since a smoother flow is formed which 
reduces the evacuation time at the expense of an increas- 
ing number of conflicts in the early stage. Furthermore 
the number of dangerous conflicts close to the door — 
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whic h hav e a direct influence on T (see the discussion in 
Sec. [VE ) — is reduced. 

Comparing the cases 2 and 3, the existence of the col- 
umn reduces the evacuation time in both cases due to the 
above reason. However, for case 2 of a central column 
evacuation times are larger than for the shifted one of 
case 3. In case 2, all the people in the room have to avoid 
the column and take a detour, while for a slightly shifted 
column some still can reach the exit directly. Thus a cen- 
tral column has a stronger screening effect which leads 
to the increase of the evacuation time compared to the 
slightly shifted case. In case 4, the column does not con- 
tribute so much to the subdivision of the flow, but acts 
more like an obstacle. For larger shifts the results will 
approach those of case 1. 



V. CONCLUSIONS 

We have extended a recently introduced stochastic cel- 
lular automaton for pedestrian dynamics by incorporat- 
ing friction effects. These are closely related to the oc- 
curance of conflicts, i.e. situations where several people 
try to occupy the same space. Such conflicts are only 
resolved with probability 1 — /i, where /i is called fric- 
tion parameter. This implies that friction as introduced 
here does not influence the motion of a single pedestrian 
(e.g. by reducing her average speed), but appears only 
through interactions. Friction effects become important 
in situations where locally high density regions occur. It 
therefore can have a strong impact on global quantities 
like evacuation times although it acts only locally. 

To elucidate friction effects and the role of conflicts 
we have investigated a simple evacuation scenario, espe- 
cially the evacuation time and the time evolution of the 
process. Without friction three different regimes with dif- 
ferent behaviour can be distinguished, as shown in [pT| . 
Introducing /i > affects these regimes in a different way. 

In general, the ordered regime, characterized by strong 
coupling to the static field and weak coupling to the 
dynamic field, is affected most. This is not surprising 
since for sufficiently large initial densities large queues 
are formed at exits and clogging becomes relevant. The 
introduction of friction leads to a strong increase of the 
evacuation time (fig. ^ due to a strong increase of con- 
flicts (flgs. |l^, ^0|). Here it is important that the con- 
flicts occur mainly very close to the exit and thus have 
an immediate influence on the outflow. In contrast, in 
the other regimes conflicts typically occur everywhere in 
the system and have therefore a weaker influence on the 
outflow properties. 

Apart form evacuation times we have also investigated 
the time evolution of an evacuation process. After the in- 
troduction of friction the outflow is irregular with periods 
of stucking. This leads to plateaus in the number {N{t)) 
of evacuated persons after time t very similar to the be- 
haviour observed in granular flow 0, |[ |9[ |l^, |ll|, |l^, ^ . 
Here the formation and breaking of arches is responsi- 



ble for the irregular flow. Similar in our case friction is 
responsible for the formation of arch-like structures that 
block the flow close to the door. The plateaus occur ran- 
domly and can therefore not be seen after averaging over 
different samples. The averaged iV(t)-curves typically 
show a linear behaviour. Using simple combinatorial ar- 
guments we found N{t) ~ • t for the /x-dependence 
of the asymptotic behaviour of the number of evacutated 
persons. This expression agrees quite well with the nu- 
merical data for densities not too small up to intermedi- 
ate values of /i. 



Finally we have shown in Sec. IVF that the introduc- 
tion of friction is essential to reproduce the surprising 
effects observed in Q when placing an additional col- 
umn in front of the exit. As in |^ this column does not 
necessarily act as an obstacle, but can — under certain 
conditions, e.g. a slightly asymmetric position — improve 
evacuation times considerably. The fact that the model 
used here differs in many respects from the social-force 
model |2l| (e.g. sign and type of interactions) used in ||] 
implies a certain robustness of this phenomenon. 
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APPENDIX A 



Defining the function 

[t/2] 
n=0 



t — n 
n 



(Al) 



we have 



[t/2] 
n=l 

[t/2] . 
ri=0 



n 



t-n-1 
n — 1 



[t/2] 



n=l 



71-1 



z". (A2) 



The first term of the r.h.s in (A2) can be rewritten as 

[(t-l)/2] 



E 

ji=0 



t-n-1 



z" = a(t - 1). (A3) 



For odd t this is obvious since then [(f— 1)/2] = [t/2]. For 
even t we have [{t — l)/2] = [t/2] — 1, but the additional 

term in the sum vanishes. 

The second term of the r.h.s in ([A^) becomes 



[(t-2)/2] 

^ E 



t-n-2 



(A4) 
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Thus we obtain the recursive relation 

a{t) = a{t- 1) + za{t-2). 



(A5) 



The solution of this recursion for a(0) = a(l) = 1 is given 
by 



a(<) 



1 



(A6) 



where x± are the roots of the quadratic equation — 
X — z — 0, i.e. 



x± 



1± Vl + 4z 



(A7) 



Therefore we have 



lt/2] 



n=0 



i + VTTliV /i-\/TT4i 



t+11 



(A8) 
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